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The critical line of the binary mixture methane-ethane is calculated via the extended corresponding states 
Van der Waals one fluid theory. The Cibbs free energy criticality criteria are solved numerically. The numerical 
derivatives are compared with the exact analytical results derived previously for the special case of the shape fac- 
tors of the extended corresponding states set equal to unity. Binary interaction parameters are adjusted to give a 
best fit of the critical line to experimental data. These interaction parameters are then used to evaluate vapor 
liquid equilibrium data away from the critical region. It appears that a fit of the critical line is not sufficient to 
obtain binary interaction parameters of general applicability. Optimization of the critical point predictions for 
the pure components is also discussed. 

Key words: Binary interaction parameters; criticality criteria; extended corresponding states; gas-liquid critical 
line; one fluid theory; Van der Waals theory; VLE prediction. 



1. Introduction feedstocks; phase equilibria is a major factor. But it is well 

known that the prediction, even the correlation, of the prop- 
The prediction of phase equilibria is both a classical erties of the appropriate systems can be exceptionally dif- 
problem of the theory of liquids and a problem of engineer- f lcu h if the results are required to any reasonable accuracy, 
ing concern. Today the chemical and fuel industries have to Prediction techniques are needed especially because the 
increase productivity and conservation and transfer to new number of possible systems makes measurement an over- 
whelming task. Prediction requires an understanding of 
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be identified clearly. This method is extended correspond- 
ing states. Here we apply it to a system of methane and 
ethane. A specific objective is to calculate the gas/liquid 
critical line and to observe the effect of the binary interac- 
tion parameters on the calculation. It is then interesting to 
see how these parameters, optimized for the critical line, 
represent vapor liquid equilibrium (VLE) data. 

The critical line in a binary mixture may be calculated by 
solving the equations 



T = T/f aa/> and V = V/h D 



(1) 
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for a temperature (T) and pressure (p) with the mole frac- 
tion (x) specified. G is the molar configurational Cibbs free 
energy of the mixture. In this work these second and third 
order derivatives were evaluated numerically, but have been 
compared with the analytical results of Wielopolski [l] 1 in 
the special case when the extended corresponding states 
shape factors are unity. The accuracy of the approach has 
thus been evaluated. 

The system methane/ethane was selected for comparison 
with experiment since the VLE data have been evaluated 
for thermodynamic consistency by Hiza, et al. [2], The pro- 
cedure is quite general, however, and we have applied it to 
several mixtures. Variations have been reported extensively 
by Watson and Rowlinson [3], Cunning and Rowlinson [4], 
Teja and Rowlinson [5], Mollerup and Rowlinson [6], and 
Mollerup [7, 8]. The overall objective is to develop a general 
technique for calculating the critical line of a binary mix- 
ture and to see if the binary interaction parameters can be 
reliably evaluated by adjusting them to give the best least 
squares fit of the critical line data. 

2. Corresponding states and equations 

The basic postulate of the theory used here — the Van der 
Waals one fluid theory — is that if the components a (a = 
l,n where n is the total number of species) of a mixture 
separately obey classical corresponding states, then their 
mixture will also obey corresponding states as if it were a 
single substance. The components can be represented by 
selected parameters, e.g., critical temperature (T$ and crit- 
ical molar volume {Vi), and the hypothetical equivalent 
substance, designed by subscript x, can be characterized by 
some suitable composition dependent averaged parameters 
T% and V%. The method then assumes that the properties of 
a pure substance at p and T, or V and T, can be evaluated 
with respect to those of a reference fluid, designated by 
subscript o, via 



where the scaling ratios h and /are defined respectively by 

/W = V'J V% and/ a „ = Tl/Tl (2) 

For a mixture the most natural definition of h z and/, follows 
from the work of Henderson and Leonard [9] to give the Van 
der Waals one fluid mixing rules: 



^z^fzfi — ^^xjt f a 0, o h a0ja 
h Xfi = El^XaXp h a8r0 



(3) 
(4) 



The cross coefficients/,^ and h a8o are left unspecified until 
further combination rules are defined, e.g., 

faPfi = Serf* (/.OMj/ftW'* W) 

h^ = Va0 , o [fc (h aa , a )» + V2 (V.)*] 3 (6) 

where fc^ and t;^ are the binary interaction coefficients 
which, although formally close to unity, can play a major 
role in the calculation of phase equilibria. 

One Fluid Mixture Equations: 

The properties of a mixture can be evaluated in terms of 
the reference substance and the ratios of eq (2). The basic 
equations are: 

Compressibility factor, Z 

Z(T,V,x) = Z£T,V,x) (7) 

= Z.IT/U, V/hJ (8) 
Molar configurational Helmholtz free energy, A 

A ( V, T,x) = A x ( V, T,x) + RTLx a lnx a (9) 
where 

A x {VJ,x) =f„ A (V/h„, T/f x J - RTtn h„ (10) 
or the molar configurational Gibbs free energy, G 

G (p,T,x) = G x (p f T,x) + RTX Xa In x a (1 1) 



where 



1 Numbers in brackets refer to the literature references at the end of this paper. 



G x (p,T,x) =/„ G (phJf XJ0 , T/f^) - RTin h xfi (12) 
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The symbol ~ refers to the molar quantity. Equations 
(7) — (12) which define the properties of an n-component 
mixture, can also be used for pure component properties if 
all subscript x's are replaced with a's. 

VLE Equations: 

For pure component VLE, equating the molar Gibbs free 
energy of each phase results in the following expression: 

[l R "IRT - (n Z - 1 + z\ ap = 

[a r "IRT - in Z„ - 1 + l\ lq (13) 

AT m A(T/f„„ W/O - A pg (T/f aa ^, V/h aa J (14) 

In eq (13) superscript_/?es refers to the residual value de- 
fined by eq (14) with A M the value of the equivalent perfect 
gas. Equation (13) is expressed in terms of the residual 
Helmholtz free energy rather than Gibbs since the reference 
equation of state has Tand F(not T and p) as the independ- 
ent variables. 

For mixture VLE one can calculate the K-value for, say, 
species a at T and p: 



(d 3 (GJRT)ldx$ T , p + 



Kc = yjx a 



(15) 



where one can derive 

RT tn K a = 0*2*%,, - fcrW (16) 

with n R " the residual chemical potential. Further manipula- 
tions give fi R " in terms of G and, for a binary mixture, 



&> = C, - J?r^ (^^\ -RTtn^ 



where 



- (nZ e -l + Z a + In 



kT 



kT 
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(18) 



RT RT B 
Critical Criteria: 
The conditions for a critical point at T,p for a mixture are 

O'SaxTr, = (d*Gidx% p = o. (19) 

Substitution of the one fluid equations gives 



(d 2 (GJRTydxX, 



1 



Xa%0 



= 



(20) 



{xjcef 



= 



(21) 



which can thus be evaluated using eq (18). 

The above equations and others have been discussed in 
full and derived by Rowlinson and Watson [3], by Eaton [10] 
and by other authors so it has been sufficient to be very 
brief. The equations form the basics of the evaluation of 
phase equilibria for a pure fluid or mixture, given the 
reference equation of state and the reference G a or A . 

Extended Corresponding States: 

In general, since classical corresponding states is not 
obeyed, equations (8) and (10) or (8) and (12) are not 
satisfied with the scaling ratios of eq (2). It is possible, 
however, to define a corresponding states so that eqs (8) and 
(10) are satisfied exactly. To do this we define shape factors 
6 and 4> so that (for a pure, for example) 



= F(~-) 



1 



n' e a 



v = r(-g-) 



hence the ratios /and h become 

Jaaja \ 7p^ ) "aap 5 " 



Vi'4>. 



= (-£-)*. 



(22) 



(23) 



The point about this redefinition, i.e., the basis of extended 
corresponding states theory, is that the corresponding 
states equations can be used formally with the provision 
that the scaling ratios are given by eq (23). It should be 
stressed that the ratios could be solved for either a pure or a 
mixture via eqs (8) and (10) but to do this would require a 
complete description of the fluids in question: essentially an 
impossibility. It is convenient to have some generalized 
analytical relation for 8 and <j>. Leach and Leland proposed 
the following [11]: 

Bo.* (Tl, K, oO = 1 + («. - u ) F[TZ, K) (24) 



0-^(72, VI, «„) = (1 + («« - u.) GiTZ, VZ)\ 



(25) 



and 



where 

TO, K) = a l + bJnTl+ (c, + dJTZWS ~ 0.5) (26) 

and 

G(r a , K) = a 2 (K + b 2 ) + c 2 {K + d 2 )in Tl (27) 

Here a> is the pitzer acentric factor or some chosen parame- 
ter and a, b, c, d are constants: 
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0.0892, 


a 2 = 0.3903 


-0.8493, 


b 2 = -1.0177 


0.3063, 


c 2 = -0.9462 


-0.4506, 


d 2 = -0.7663 



a y = 

*i = 

c i = 

d x = 



The asterisk denotes the value reduced by the critical value. 
The equations are constrained in that V% is set equal to 2.0 
for all VI > 2.0 and to 0.5 for V* a < 0.5: T£ is set to 2.0 if 
Tt > 2.0. 

We [12] have recently tested the Leach-Leland equations 
for the hydrocarbons C,-^,, over an extensive range of ex- 
perimental conditions and revised coefficients are reported 
in the reference. We also verified that the original equa- 
tions were satisfactory for reduced temperatures greater 
than 0.5. 

3. Calculation and numerical methods 

The objective is to solve the critical criteria eqs (20) and 
(21) for the methane/ethane system and in so doing, observe 
the effects of the interaction parameters £ and tj of eqs (5) 
and (6) on the results. Having these values, we then evaluate 
some -K-values for selected temperatures using eqs (15)— (18). 
We chose methane as the reference fluid, the equation of 
state for which is the 32 term BWR of McCarty [13]. Critical 
parameters and Leach-Leland acentric factors for methane 
and ethane are given in table 1. 

TABLE 1. Parameters for Methane and Ethane 



where 



7* 
(K) 



(cm'/mole) 



P' 
(Bar) 



CH 4 



190.555 
305.33 



97.75 
147.06 



44.793 
47.448 



.005 
.105 



Analytical and Numerical Evaluation of the Derivatives: 

The numerical techniques used in this work are standard. 
We use the central difference formulas [14] for which the 
first two terms in the infinite power series expansions are 
given here. For the derivatives of a function/evaluated at a 
point x, one has 



4A 

dx), 



2h \2h 

*f\ .. /.-2/ +A / 2 -4/ 1 -f6/ -4/.,-f/- 2 



dx* ),~ 



d 3 f \ _ 

a* 3 A ~ 



\2h 2 

/ 2 -2/ 1 -f2/- 1 -/- 3 
2h 3 8h 3 



/ 3 -4/;+5/;-5A+4./W- s 

8h 3 



(28) 



(29) 



(30) 



f.=M 



f n = f{x + nh) 



(3D 



The difficulty is to choose a value of h which is not too small 
(otherwise significant figures will be lost in evaluating the 
numerators of eqs (28)-(30)) but not too large (otherwise the 
truncation error, which can be estimated by the second term 
of eqs (28)-(30), will be large). One also has to consider the 
word length of the computer and the convenience of using 
single versus double precision. In this work we calculated on 
a CDC 6400 and a CDC 6600 machine with a 60 bit word 
length (13 significant figures). 

We were able to observe definitely the effect of varying h 
for the special case = = 1, i.e., for classical correspond- 
ing states. Equations (19) and (20) have been solved analyt- 
ically by Wielopolski (1980) and the lengthy expressions are 
reported in an NBS publication [1] and will not be repeated 
here. For example, table 2 lists the number of figures in the 
numerical results which were in agreement with the analyti- 
cal results for the first, second, and third derivatives of 
GJRT for a particular test case. The number of figures in 
agreement for the function value of GJRT itself was 10-12. 

TABLE 2. Comparison of Numerical and Analytical Results for Derivative 

Calculations Using Single Precision Arithmetic. Subscripts refer to the 

derivation with respect to x. 



(GJRT), 



{GJRTh, 



{G,IRT) 3r 



10" 
10" 

io- 



Table 2 indicates that the first order derivative is trunca- 
tion error controlled, since its value becomes more accurate 
as h is decreased. The second and third derivatives are, on 
the other hand, controlled by the loss of significant figures 
since as h is decreased, they lose accuracy. Since the third 
order derivative is the least accurate, we chose the value of 
h for which it is calculated most accurately. 

We now consider what the smallest values of the second 
and third order derivatives are which can be calculated with 
h = 10" 3 , since our eventual goal is to solve the equations 
for the critical point by driving the values of those 
derivatives to zero. The derivatives go to zero by a cancella- 
tion of the two terms in eqs (20) and (21), that is, the con- 
tribution from the hypothetical substance is cancelled by 
the ideal mixture contribution. For this reason, the values of 
the derivatives cannot be made arbitrarily small. The ideal 
mixture contribution (which can be computed with negligi- 
ble error) can only cancel as many significant figures as ap- 
pear in the hypothetical substance contribution. Consider 
the case in table 2 with h = 10" 3 . For the second order 
derivative, the hypothetical substance contribution has five 
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significant figures, and its value is order unity (abbreviated 
0(1)). If the ideal mixture contribution were to cancel all five 
of these figures, the result would be a number of 0(10" 5 ) with 
no significant figures remaining. For the third order deriva- 
tive, the hypothetical substance contribution contains four 
significant figures, and is 0(10). Cancelling all significant 
figures would leave a number of 0(10~ 3 ). 

In our first attempt at calculating critical lines based on 
the numerical evaluation of the derivatives in eqs (20) and 
(21) using single precision arithmetic, we were unable to ob- 
tain convergence of the temperature and pressure to five 
significant figures. The problem appeared to be that there 
were not enough significant figures in the derivative calcu- 
lations. While the truncation error is inherent to the for- 
mulas being used, the loss of significant figures can be com- 
pensated by adding more figures to the function values. 
This was done by the use of double precision arithmetic 
which gives us 26 significant figures on the CDC 6400. In 
table 3 below, the results for the numerical derivatives 
calculated using double precision arithmetic are compared 
with the results arrived at analytically. Again, reported in 
the table are the number of figures of agreement between 
the two results. 

TABLE 3. Comparison of Numerical and Analytical Results for Derivative 
Calculations Using Double Precision Arithmetic. 



(GJRT), 



(GJRTk, 



(G,IRT) 3x 



io- 

10" 
10" 
10" 
10" 
10- 
10" 
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5 
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9 


6 


7 


8 
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7 


8 


7 


7 


8 


7 


7 


8 


4 


7 


8 


1 



For h = 10" 3 , the single and double precision results are 
the same, which indicates that truncation error is control- 
ling. Looking at the double precision results, the third 
derivative shows an increase in accuracy as h is decreased to 
10" 5 ; clearly indicating that the truncation error is decreas- 
ing to this point. As h is decreased past 10" 7 , accuracy is lost 
due to loss of significant figures. 

Based on these results, a value of h = 10" 6 is chosen to 
compute the derivatives in double precision. Given this val- 
ue for h, the smallest value of the second derivative which 
may be calculated (containing no significant figures) is 
0(10- 8 ), and that for the third derivative is OflO -6 ). The 
calculations of the binary critical line were subsequently 
made to converge to five significant figures for both temper- 
ature and pressure. 

4. Results 

It must again be stressed that the general procedure for 



calculating the critical line or VLE is predictive and re- 
quires only the critical constants and an acentric factor for 
the fluid of interest, or of the components in a mixture. For 
a relatively simple system the results will be reasonable 
without optimization of any parameters. Since, however, we 
are concerned only with VLE and the critical point we con- 
sidered two straightforward optimization procedures involv- 
ing the factor w. The first was to adjust w to give the best 
representation of the pure component vapor pressure curve, 
the second was to force the critical temperature and 
pressure of the pure fluids to correspond exactly with those 
of the reference substance. This second variation is simply 
to set w a = w : hence by eqs (24) and (25) = 1 and <f> = 
Zl/Za—a form of classical corresponding states. One should 
note that the two procedures are not the same because the 
Leach-Leland equations are not constrained at the critical 
point. 

Ethane: Pure Component Results: 

We first considered the ethane vapor pressure curve 
which was obtained using Leach's expression for the shape 
factors. The value of the acentric factor for the Leach equa- 
tions was determined by optimizing agreement with the 
vapor pressure data by a trial and error procedure in which 
the sum of the average absolute deviations, for the vapor 
pressure, and saturated vapor and liquid densities, were 
minimized; temperature being chosen as the independent 
variable. The temperature range over which the results were 
optimized was 180 K to 300 K; 180 K being the lower limit 
for which Leach's equations were designed; 300 K corre- 
sponded to the maximum temperature for which the vapor 
pressure program would converge. Calculated and data 
values were compared at 10 K increments. The average 
deviations obtained for several values of acentric factor are 
given in table 4. The value chosen for acentric factor in this 
work (w = 0.094) is seen to give a substantial improvement 
over the Pitzer value (w = 0.105) which was used by Leach. 



TABLE 4. Variation of Ethane Vapor Pressure Curve With Acentric Factor 
(Leach 6, <j>). 



Acentric 
Factor u 


Ave % 


Ave % 


Ave % 




|Ap| 
P 


|Ae~*l 


|Ae«.l 

Qii, 


E Ave % 


0.105 


2.080 


3.148 


0.6158 


5.844 


.100 


1.259 


2.215 


.3706 


3.845 


.096 


0.595 


1.460 


.3791 


2.434 


.095 


.467 


1.270 


.3894 


2.126 


.094 


.444 


1.164 


.4067 


2.015 


.093 


.479 


1.127 


.4245 


2.031 


.092 


.543 


1.137 


.4424 


2.122 


.090 


.736 


1.241 


.4904 


2.467 


.005 


60.579 


61.633 


4.746 


126.958 
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The curves for vapor pressure and orthobaric densities 
(q = II V), obtained using Leach's shape factors with the 
optimized acentric factor are compared with the correla* 
tions of Goodwin, et al. [15] to obtain the deviations platted 
m figure 1. 

Also in figure I, deviations for the vapor pressure and 
orthobaric density curves predicted using classical corre- 
sponding states are presented. This figure emphasizes that 
the Leach shape factors make a significant difference. The 
vapor pressure deviations are positive, and become larger 
as the triple point is approached, since classical correspond- 
ing states predicts a slope of the vapor pressure curve 
(dP/dT) which is too small. The deviation of pressure goes 
to zero at the critical point because the two parameters are 
choosen to make the critical temperature and pressure cor- 
respond exactly. 



The critical point results for ethane are in table 5. Notice 
that the results are better using classical corresponding 
states than with the Leach shape factors. This is because 
classical corresponding states forces either the critical 
temperature and density to correspond ($ = (f> = 1), or it 
forces the critical temperature and pressure to correspond 
{8 = 1, tf, = Z'JZ% 



TaBLE 5- Eihane Critical Point Prediction* 
6=$ = \ 



Data 



Leach 0, 
#«? = 0.09-1) 



3 = 1. 



7; (K) 305,33 


307,01 (O.SS)' 


305.33 


305.33 


P, {bat) 47.488 


46,790 (2.S3) 


47.7S0 (0.SS) 


47.448 


q c (mol/L) 6.80 


6,98(2.65) 


6.80 


6.76 (-.59) 



* Percent deviation is in parentheses. 
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FiGUHE I. Petceni deviations [(e*pi-ca!c) X 100/calc] for the vapor 
pressure, vapor density and liquid density for ethane at saturation. Data 
from Good win [15]. Calculations from the e s ! ended corresponding states 
method with the parameter w„ optimized, circles. Also shown as the dashed 
line are the results with 8 = 1 and <£ = ZJ/Z* , see equations (22)— (27). 



The Critical Line: 

We first calculate the critical line using the Leach shape 
factor equations with acentric factors of 0.005 and 0.094 for 
methane and ethane respectively. The results are plotted 
against the critical line data found in the review article of 
Hicks and Young [16], and identified in the caption to 
figure 2 (the symbols used in figure 2 are identical with 
those used in figures 3 through 7). 

The results are presented in the form of T-x and p-x plots 
in figures 2-5 and show the general trends obtained by 
varying the binary interaction parameters, £ and ij. Holding 
£ constant, figures 2 and 3 show that tj has a small effect on 
the Tx curve, and a large effect on the p-x curve. In both 
cases, increasing ij gives a better representation of the data. 



€=1.00, fj=1.0Q 




220- 



o.a 



X CH 4 

FIGURE 2. The methane-ethane T-x critical plot calculated via ei tended 
corresponding slates showing the variation caused by the interaction 
parameter, n. Dala [16[: v" Bloomer, Cani and Parent; D Price and 
Kobayashi; O Wichelerje and Kobayashi.The values for the pures, A, are 
from Goodwin [15] and [19J. 
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Holding 7j constant, figures 4 and 5 show that £ has a much 
larger effect on the T-x curve than did tj, and an equally 
large effect on the p-x curve. The important point to notice 
is that the maximum value in the p-x curve is shifted 
towards small mole fraction values (of CH 4 ) by decreasing £. 
The best representation of the P-x curve in figure 3 (i.e., 
£ = 1.00, 7j = 1.08) indicates that the peak of the curve 
needs to be shifted towards the smaller mole fractions to im- 
prove the agreement, thus, £ should be decreased. 



£ = i.oo, 17=0.95 
£=1.00, 17=1.00 




*CH, 



FIGURE 3. The p-plot corresponding to figure 2. 




l CH, 



FIGURE 4. Variation of the parameter, £, for methane-ethane in the T-x 
plot. 



To achieve the goal of obtaining the interaction parame- 
ters by fitting the critical line data, a manual search tech- 
nique was initiated. The "best fit" was defined in the least 
squares sense. The results of this search were that £ = 0.97, 
and 7j = 1.13 were chosen as the "best" values for the in- 
teraction parameters. The "best fit" T-x and p-x curves are 
presented in figures 6 and 7 respectively. 



(=0.95,17=1.00 

70h 




v ch„ 



FIGURE 5. The p-x plot corresponding to figure 4. 



(=0.97,17 = 1.07 




FIGURE 6. The T-x critical line of methane-ethane with optimized interac- 
tion parameters. The two curves correspond to extended corresponding 
states [£ = 0.97, ij = 1.13] and the classical corresponding states with a)„ 
= «.K = 0.91,7} = 1.07]. 




^CH* 



FIGURE 7. The p-x plot corresponding to figure 6. 
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The fit of the T-x curve is good, with only one data point 
which seems astray. The p-x curve, however, does not have 
the right shape to fit the data well. Part of the fitting prob- 
lem is due to the bad prediction which is made for the criti- 
cal point of pure ethane. This led us to try the second ap- 
proach of setting oj = oj . Hence, the critical endpoints in 
the T-x, and p-x curves are exact. A new optimization led to 
the parameter values £ = 0.97, rj = 1.07. While the fit of 
the T-x curve was not significantly improved that for the p-x 
curve was. These results are shown in figures 6 and 7. 



Vapor-Liquid Equilibria Results: 



one is subcritical (144.26 K). The sources of the data are: 
250 K isotherm, Davalos, et al. [17]; 199.92 K and 144.26 K 
isotherms, Wichterle and Kobayashi [18]. 

The VLE calculations used the Leach shape factors with 
the acentric factors 0.005 and 0.094 for methane and ethane 
respectively. The results are presented as K-\alue deviation 
plots for both the methane and the ethane K-\a\ue predic- 
tions. Figures 8, 9, and 10 contain these curves with the in- 
teraction parameters obtained from the critical line fit (i.e., 
£ = 0.97, rj = 1.07). These figures also show that setting the 
interaction parameters to unity gives much better VLE pre- 
dictions than do the parameters obtained from the best fit 
of the critical line data. 



Of the VLE data judged to be thermodynamically consist- 
ent by Hiza, et al. [2], three representative isotherms were 
chosen to test the predictions made using the binary inter- 
action parameters determined in the previous section. Two 
of the isotherms are supercritical (250 K and 199.92 K), and 
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FIGURE 8. Methane-ethane K-value deviation plots at 250 K. Shown are the 
(dashed) curves with the interaction parameters from the critical line fit 
[{ = 0.97, 17 = 1.07] and, for reference, with £ = 17 = 1.00 as the solid 
curves. 
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FIGURE 9. Melhane-e thane iC-values at 199.92 K. 



5. Summary and Conclusions 

The proposed technique of calculating binary critical 
lines by numerically evaluating the second and third order 
derivatives of the Gibbs free energy has been checked with 
an analytical solution for the special case of classical Cor- 
responding states, and has proven successful. The best least 
squares fit of the critical line data of the system methane- 
ethane was then shown to be poor (particularly the p-x 
curve) if the Leach shape factors are used with an acentric 
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factor optimized for pure component vapor pressure predic- 
tions. This is due to a bad prediction of the critical endpoint 
for ethane. To improve this fit, we use classical correspond- 
ing states to force correspondence of the temperature and 
pressure at the critical line endpoints. However we also 
show that the pure component vapor pressure predictions 
are not satisfactory if this is done. Finally, VLE predictions 
are made using Leach shape factors with the acentric factor 
optimized for vapor pressure predictions, and the binary 
interaction parameters obtained from the best fit of the crit- 
ical line data (i.e., with 6 = 1, = Z%IZ% The results are 
not as good in general as those which are obtained by set- 
ting £ = i] = 1. Hence we conclude that a fit of the binary 
critical line does not yield binary interaction parameters of 
any general significance. 
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FIGURE 10. Methane-ethane lvalues at 144.26 K. 
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